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Abstract 

We solve the Bethe-Salpeter equation for hydrogenic bound states 
by choosing an appropriate interaction kernel K c . We want to use our 
solution to calculate up to a higher order the hydrogen Lamb-shift, 
and as a first application we present up to order (a/vr) (Za) 7 the 
contribution of the lowest order self-energy graph, calculated exactly. 
The basic formalism is a natural extension to the hydrogenic bound 
states of the one previously presented by R. Barbieri and E. Remiddi 
and used in the case of positronium. 
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1 Introduction 



Bound state systems like positronium and muonium are a good test of QED. 
Hydrogenic atoms are not completely reducible to a QED problem because 
the structure of the proton. In particular the finite size of the proton gives 
rise to a lower limit of the theoretical calculus precision (for a good review 
see [0] and more recently 0). 

However, there are two kind of considerations which make hydrogenic 
atoms still interesting from a QED point of view. First, the requested pre- 
cision is not yet achieved in many theoretical predictions, e.g. the hydrogen 
Lamb-shift requires a theoretical precision of about 1 kHz (the present sta- 
tus of the calculation of the hydrogen Lamb-shift can be found in [|]), and 
new but incomplete contributions have been recently calculated (e.g. ||, @). 
Second, the fully detailed treatment of a relatively simpler problem as the 
hydrogenic atom (which is reducible in QED to a one particle problem in 
an external potential) can give a hint for the treatment of more complicated 
system as the positronium. For this purpose it is necessary to use the same 
formalism to describe the positronium and the hydrogenic atoms. 

In sect. 2 we extend to the hydrogenic atoms the formalism which was 
proposed by R. Barbieri and E. Remiddi M for the positronium (which we 
call BR formalism). Following || we propose as interaction kernel a sort of 
relativistic dressed Coulomb interaction, so that the Bethe-Salpeter equation, 
which for bare Coulomb interaction is the Dirac-Coulomb equation, is ana- 
lytically solvable in closed form. Then we write a perturbative expansion for 
the energy levels which immediately reproduces the Dirac levels. In sect. 3, 
as a first application of our formalism, we calculate analytically up to order 
a/ir(Za) 7 for the levels n = 1,2 the first contribution to the self-energy to- 
gether with other graphs necessary to cancel spurious terms, which typically 
arise in this sort of calculation. An appendix is devoted to review and discuss 
the method proposed by J. Sapirstein [0] to treat perturbatively the Dirac 
propagator. 
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2 BR formalism in the context of the hydro- 
genic atoms 

In QED hydrogenic atoms are well approximated by an electron moving in 
an external field V c (r) = —Za/r, that represents his Coulomb interaction 
with the nucleus. The Green's function G(W;p,q) of this electron contains 
all the necessary informations about the system, and it can be written in 
perturbation theory as the sum of all the Feynman's graphs with an incoming 
and an outcoming electron leg (see fig. 1; W is the energy and p,q are the 
incoming and outcoming spatial components of the electron's momentum; 
the ball includes all radiative corrections, the external lines are Coulomb 
interaction vertices). 



Fig.l The Green's function G(W; p, q) of an electron in external field. 

It is well known that in bound state theory the Green function G has 
simple poles in the energy W. For hydrogenic atoms these poles are grouped 
round the Balmer levels —m(Za) 2 /2n 2 and are labeled by the quantum num- 
bers n, I (the angular momentum which is an " almost good quantum num- 
ber" in the sense of |§) and j (the total momentum). For every pole W n ij 
we can write: 



(W,p ) 




(W,q ) 



X 



X 



G(W;p,q) 



Rnij(p,q) 
w - w nlj 
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where R n ij is the residuum and G n ij is the regular part of G at W = W n ij- 

The standard way to obtain a perturbative expansion of W n ij starts from 
the Bethe-Salpeter equation M: 



G(W;p,q) 



+ 



G (W;p) 
_fk_ 

(27T) 5 



(27i) 3 5 i3 \p-q) 
K(W;p } k)G(W;k,q) 



(2.2) 



where Gq(W;p) = i(^ Q W —p- ; y — m+ie)~ 1 is the electron free propagator and 
K(W;p,q) is the interaction kernel (i.e. the sum of all 1-particle irreducible 
graphs with external fermionic legs removed). Then we write: 



K{W- p, q) = K C (W; p, q) + SK(W; p, q) , 



(2.3) 



which is, rigorously speaking, a definition of 5K, once K c is suitably chosen. 
To choose K c , we ask that it satisfies the two conditions: 

(i) in the non relativistic limit K c — > —ijo Vc, 

(ii) the Bethe-Salpeter equation for K c : 



G c (W;p,q) 



+ 



G (W;p) 



{27i) 3 5 i3) {p-q) 
K c (W;p,k)G c (W;k,q) 



(2.4) 



can be explicitly solved in closed form. 

From condition (i) it follows that G c can be considered as a zeroth order 
approximation of G. Therefore G c has singularities in W = W£ [] and W£ ~ 

m — m(Za) 2 /2n 2 : 



G c (W; P ,q)= > 



G c n (W;p,q), 



(2.5) 



where R c n ij is the residuum and G c n is the regular part of G c at W = 
From condition (ii) it follows that W£, R^ij an d G c n are explicitly known. 

1 According to general expectations the exact problem has less degeneracy than the 
unperturbed one. 
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Using the quantities denned above, first we can write the Bethe-Salpeter 
equation ( |2.2|) in terms of G c and SK as: 



G(W;p,q) 



+ 



G c (W;p) 
d 3 k 



SK(W;p,k)G(W;k } q) 



(2.6) 



and then the perturbative expansion of the energy levels: 

W nl] = W: + ±Tr{5K(W:)R c nlj } 

+ ^Tr{5K(W:)G c n (W:)5K(W:)R c nl3 } 



(2.7) 



where we have omitted for simplicity the explicit indication of the momenta, 
and D is defined as the following trace 0: 



D = Tr 







:G-\W c n )R c ri 



(2.8) 



K dW~ c y " n ' nl if' 

Each term of ( |2.7|) is a series in (Za) with fixed (a /it). We observe that 
the expansion fl2.7|) is "perturbative" in the sense of increasing orders of 
8K(W£). For consistency the explicit calculation must exhibit that to an 
increasing order in 8K(W£) it corresponds an increasing leading order in 
Za. In ]TD| the reader can find a more detailed discussion. 
In |[L1| , (see also it was found a kernel if c satisfying conditions (i) 



From (^J) it follows that G c = G Q — K c and then we can write ( |2.8| ) also as: 
£> = -Tr {i^Rlu } - Tr ^K c (W^)R c nlj \ . 
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and (ii) for positronium. In what follows we choose for K c []: 



K c (W;p,q) = -i f(W;p,q)V c (p-^°A + (f) 



1 + 7° 



(2.9) 



where V c \ 
and 



A±(p) 



4nZa/p 2 is the Fourier transform of the Coulomb potential, 
E p ± (m — p ■ 7)7° 



2 -En 



f(W;p,q) 



En 



16m 2 E n E, 



(Ep + m)(E p + W)(E q + m)(E q + W) j 



p 2 + m 2 



In the non relativistic limit p, q — > it is easily seen that K c satisfies 
condition (i). 

If we define H c (W;p,q) as: 



G c {W;p,q) = G (W;p) 



(27r) 3 5 (3) (p-g) 



z/(W;p, g )7 A + (p1^^A + ®G (W;g)F c (W;p,g) 



(2.10) 



from the comparison of ( |2.1(J| ) and ( ^.4[ ) we obtain the following equation for 



H c (W;p,q) = V c (p-q) 
d 3 k 

+ 



(2.11) 



(27r) 3 fc^_|i +ie 

2m 2m 



V c {p-k)H c {W; k,q); 



3 An other choice can be K c = Kd = —i^°V c . As a consequence of this choice G c = Qd 
(where gD is the Dirac-Coulomb propagator and will be defined later), go is known in 
analytical closed form and its poles and residuum at the poles also (see |l3| and for a 
detailed study |H|]). But analytical computation with these quantities are more complicate 
than in our formalism because the unpractical structure of go- 
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the solution of ( [2.11| ) is known and can be written in the Schwinger integral 
representation |15|: 



H c (W;p,q) = V c {p-q) + (Zaf 



Arcm 



(p — q) 2 \/m 2 - W 2 



dp 



P 



— mZa 
y/m 2 ~W 2 



(E^-W^Ef-W 2 ), _ W 
4(m 2 ~W 2 )(p-q} 2 V 1 V) 



(2.12) 



Inserting ( 2.12|) into ( p. ICQ we obtain: 
G c (W;p,q) = (27r) 3 5^(p-q)G (W;p) 
+ i 



f{W ^ q) A + (p)^A +(g -) 7 0K(p-,-) 



(E p -W)(E q -W) 
1 + Za 



m 



1 n\/m 2 -W 2 

dp- ' 



(2.13) 

so that condition (ii) is also fulfilled. 

From ( 2.13f) ) we see how the Green function G c has poles at the values of W 



W c = m\\- 



x {Za? 



rr 



(2.14) 



and isolating the singular part from the regular one at W = W£ in ( |2.13|) we 
find also the residuum R-nij a ^ the P°l e: 



W c 



E p + m 



1 +7 ! 



n-l 



Eg + m J 
21 + 1 (p- q 



1=0 



pq 



A + (<?>°, 
(2.15) 



where Pi(z) is the Legendre polynomial of order / and R n i the radial part of 
the Schrodinger- Coulomb wave functions. 
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G^(W;p,q), the regular part of G c at W = W£, can be obtained explic- 
itly subtracting from ( j2.13|) the singular part of ( |2.5|) and taking ( |2.15|) into 
account (see [IB] for the explicit case of positronium). 

To obtain from fl2.15|) the expression of R^ij, following |1(J, we write the 
identity: 



where 



21 + 1 n (p-q 



Air 



-Pi 



pq 



p q 

y 



p q 

y 



p-q 

k vq , 



+ 



j - i ( p^jq^j + Ell\ — 

2n \ p q pq I dz 



Pi(z, 



(2.16) 



R^ij is then given by: 



RC (if * * (E P (Z P + TO V (E q {E q + W< 

Rnlj{P,q) = yy- c \ 



E v + m 



E a + m 



R nl (p)R nl (q)K + ^ l -±^D{ 



p q 

y q, 



which is an eigenfunction of parity and total momentum. 

d 3 k 



R^ij satisfies: 



RUp,q) = G (WZ;p) 



(27T) 



-K c (W:;p,k)R^(k,q) 



(2.17) 



(2.18) 



known as the Bethe-Salpeter equation for the residuum. 

The expansion (|2.7|) is now completely defined. The first terms of the 
expansion that must be calculated are D and Tr ^K c {W^)R c nl ^. These terms 
are not related to any Feynman graph, they only depend on the formalism 
which we have adopted. 
After explicit calculation we find: 



D = 2j + 1, 



(2.19) 



7 



and 

- ^ {Kdw^} = ~wr^~- (2 - 20) 

D is the degeneration of the level (n, l,j) Q and (K c ) n is, at the leading order 
(apart for a factor 2), the Balmer series (which follows from condition (i) on 
K c ). 

The next term of the expansion (|2.7|) giving the leading corrections to 
the level is Tr ^KiW^jR^^. In fig. 2 we show K up to one loop radiative 
corrections. 



X 



Fig. 2 The first terms of the irreducible kernel K for hydrogenic atoms. 



In fig. 2 with the doubled line we have represented the sum of graphs of 
fig. 3. These graphs can be resummed using the Dirac-Coulomb equation: 



g D {W;p,q) 



+ 



G (W;p) 
d 3 k 



(27T) 



-iV c (p-k)j°) g D (W;k,q) 



(2.21) 



This equation is of the kind of the Bethe-Salpeter equation ( |2.2| ) with kernel: 

K D (p 1 q) = -zV c (p-q) 1 °. (2.22) 

A graphical representation of is given by the first graph of fig. 2. Therefore 
equation ( |2.21| ) can be treated perturbatively in the formalism of equations 



This result less than surprising is a natural consequence of the Bethe-Salpeter 
formalism. 
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(|2.2|)-( |2lf ) with K c given by equation ( |2.9| ) and consequently SK = Ke> — K c . 
An alternative approach to treat perturbatively the Dirac-Coulomb propa- 
gator g D was given by J. Sapirstein |7J and can be found in appendix. 



X 



X X 



Fig. 3 The Dirac-Coulomb propagator go- 



The contribution to the energy levels of K D originates up to order (Za) 4 

W D : 

1 r d 3 p d 3 q 



the well-known Dirac levels WE'a 



{Kd) 



nl i DJ (27r) 3 (27r) 3 

2 / J 



Tr{-iV c (p-q)j°R c nlj (p,<f)} 



(Za) 2 (Za 
—m m 



2n 3 \j + 



1 I > 

2 / 



ir, 



D 
"J 



W: + (K D ) nlj - (K c ) n 

(Za) 2 (Za) 4 ( 1 
m — m m 



2n 2 



2n 3 \ j + \ An) ' 



(2.23) 



(2.24) 



More interesting are the one loop corrections. In particular in the next sec- 
tion we discuss the contribution to ( ^7|) of the second graph of fig. 2, the 
self-energy graph. 



3 Application to the self-energy contribution 

The first one-loop correction to the hydrogenic atoms' energy levels is given 
by the self-energy graph (second graph of fig.2; a recent calculation of this 
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contribution in an other method can be found in fl7|| ). 

Let us recall that the contribution to the self-energy kernel due to the 
first graph of the right hand side of fig. 4 could be written as: 

ffi « = / g^ vAW -^fcfe (3 ' 1) 



X X 



Fig. 4 The (Za) expansion for the self-energy graph. 



Expression ( [3.1|) is obviously divergent. Choosing the Pauli-Villars regu- 
larization one has: 



i S reff (p; A) 



d A k 



-ig^ v -A 2 



k 2 + ie k 2 - A 2 + ie 
a , 



i + m)A(p 2 ; A) + m B(p 2 ; A)], 

7T 



(3.2) 



where 



1 



1 / m 2 



A(p 2 ;A) = -(l-m 2 /p 2 )(l + m 2 /p 2 )\og(l-p 2 /m 2 )--[ — -l 



4 

51 ( A 2 
8 4 ° g l^ ' 
1 



4 V p 2 



1 ( m 2 



B(p 2 ;A) = -(l-m 2 /p 2 )(3-m 2 /p 2 )\og(l-p 2 /m 2 ) + -[ — -l 



4 

3_3 b /A 2 ' 
8 4 ° g \ m 2 



4 V 



(3.3) 
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The mass-shell renormalization prescription is imposed by defining the mass 
renormalization counterterm: 



i 5m = — iS re9 (p; A) 



so that 13 



% 5m = i—m 



3 3, 

8 + 4^ 



(A- 

\ m? 



(3.4) 



(3.5) 



The relevant kernel for our purpose is the regularized and mass-subtracted 
graph: 

K s . e .{p; A) = z {Kegip; A) + 5m) , (3.6) 

which is still u.v. divergent (wave-function renormalization has not been 
carried out). If we rearrange the terms to isolate the A dependent part Kf~^ 
from the rest P , we will obtain: 



where 



K s . e .(p;A) = K«Z(p;A) + K s . e .(p) 



K2(p;A) 



.a 

7T 



m) 



a 



71 



G£ l (p) --tot 



- log ( A_ 

4 \m z / 

m 2 / / ' 



(3.7) 



(3.8) 



and 



K s . e Xp;A) 
A(p 2 ) 



B( P 2 ) 



+ 



« ~[(-^ + m)A(p 2 ) + m i?(p 2 )], 



P 



1 frn 2 
4 I p 2 



m 



4 I 1 p 2 



1 / m 
4 I p 2 



- 1 



m 



1 + — log(l-p 2 /m 2 ) 



P 



m 



3- — \\og(l-p 2 /m 2 ) 



P 



(3.9) 
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The contribution of ( |3.7|) to the energy levels according to (|2.7|) is given by: 



nlj \ ' I nlj 

(3.10) 



Because R c nl j satisfies the Bethe-Salpeter equation (|2.18| ) and taking in ac- 
count equation (|3.8|) we can easily calculate (Kf v 



s.e. i j ■ 




The explicit value of (K c ) n has been given in ( |2.20| ). 

We don't need to perform a wave-function renormalization: if we adopt 
a gauge invariant regularization scheme the wave-function renormalization 
counterterm is exactly compensated by the vertex renormalization countert- 
erm, hence one obtains the correct physical result without taking in ac- 
count these counterterms 5 . The A dependent term (|3.11|) must also cancel 



if we evaluate the vertex contribution to the energy expansion ( |2.7| ) in the 
same gauge invariant Pauli-Villars regularization scheme adopted for the self- 
energy and without subtracting the renormalization counterterm (see for a 
detailed discussion in the positronium context fT2|j ). 

We have therefore to consider the second graph of the right hand side of 
fig.4. The kernel due to this graph is: 

K ver = -iV c T° reg , (3.12) 

where the Pauli-Villars regularized T^ eg , 

d A k , 



r? eg (p,q;A) = J 



7^l^~ ie ^h — u ~^^~A — v —(-ieY 

(27TJ 4 f — If — m + te q — f — m + te 



° Moreover in this manner one avoids the problem of the spurious infrared divergences 
which may arise from mass-shell wave-function renormalization. 
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satisfies the Ward identity: 

(p - Q)^ eg (p, g; A) = Ke 9 (p; A) - K eg (q; A). (3.14) 



From ( p.!4j ),( |3~2D and ( |3.3|) , differentiating with respect to p u and rearranging 
terms, it follows that: 

Keg(P,<i; A ) = ^S re9 (p;A) - (p-g) M — r^(p,g;A) 



(3.15) 



where T u is u.v. finite (A independent). 

The contribution of ( |3.12|) to the energy levels is also given by: 



(Kver) nlj = ^-Tr{^V c T° reg (W:)R c nlj } (3.16) 



D 

a ( 1 , /A 2 



7T 



where we have defined K ver = — iV c T°. The explicit value of {K D ) nl - has 
been given in ( |2.23| ). 



Comparing (|3.11|) with (|3.16|) we note that only the a/7r(Z'a) 2 log(A 2 / 



m 



2\ 



terms, corresponding to the leading order in Za, cancel. Then for a complete 
cancellation of the divergent terms it is necessary to sum other divergent 
terms arising from expansion ( |2.7| ) To obtain for instance the cancellation 
up to order a/7r(Za) 4 log(A 2 /m 2 ) the terms which must be considered arise 
from the last contribution explicitly written in expansion (|2.7p: 



6 We note that this problem doesn't occur if we use the kernel K c = Kjj in this case 
the divergent part of (3.16) completely cancel ( [3.11 ). 
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Because (taking in account ( |3.8| ), ( P7T5| ) and note 2): 



one loop log 



(s) 



terms 



a 



a / 1 
- --log 

7T \ 4 



a 

7T 



-^log 



77T / 

'A 2 ; 
m 2 

m 2 



/ - K div 

\dW se -' 



I nlj 



(3.17) 



-Z7 



9 



(l + 0((Za) 2 )) 



one has: 



((Kd) 



(K D - K c ) 1 



a 



lj \dW"^\ 



5K 



nlj 



a 



A 



one loop loj 



a 



terms 



A 2 



- -7 l°g [—, +0 -(Zaf log — 



7T 



7T 



(3.18) 



Summing now (|3.11| ), (|3.16|) and (|3.18|) up to order a/n^Za) 4 all log(A 2 /m 2 ) 
dependent terms cancel. 

The contribution of ^i^ s . e .^> ^ has been evaluated on the levels n = 1, 2 

exactly, as a first step we show only the leading order terms: 



K,,.) u = ^(Zaf 



K s ,) 2 = -(Zaf 

I 2s 7T 



^4, = ^ Za ^ 
I 2p 7]- 



+ 21og(Za) + 21og(2) 



-± + \log{Za) 

11 1 

1 — log (Zen) 

96 2 Sl ; 



(3.19) 

(3.20) 
(3.21) 



Here and in the following the indication of just the nl levels means there isn't 
contribution to the splitting in the j levels. 

The a/n(Za) 2 and ct/7r(Za:) 2 log(Za!) terms in fl3.19|) -( [3.2lD are expect 
to vanish. This occur because the Dirac levels ( |2.24| ) of order (Za) 4 are 
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completely given by the graphs discussed in the previous section. This is 
a well known feature of the Feynman gauge (see e.g [0); it decreases the 
speed of convergence of the perturbative series by generating spurious terms 
of low order that only at the end of all calculation (up to requested order in 
Za) must cancel each other. 

Up to order a/n(Za) 2 the contributions to the energy levels n — 1, 2 of 



K„ 



nlj 



are: 



k ) 


a 

I s 7T 




~(Za) 

7T 


) = 


-(Za) 

7T 



iZa) 2 
1003 



39 



2\og(Za) - 101og(2) 



288 
2009 



864 



I hg{Z a) - f log(2) 

1 32 
-log(Za)--log(2) 



(3.22) 

(3.23) 
(3.24) 



These contribution only cancel the a/n(Za) 2 log(Za) spurious terms of 
( |3.19|) -( |3.21| ). To obtain the full cancellation of all aj'niZa) 2 terms which are 
in ( |3.19|) -( |3.21|) it is necessary to calculate the leading order of the remaining 
graphs of fig.4: 



ladder, 



nlj 



D 



E 

i=i 
d 3 Xi 



d 3 p f d 3 q 



(2vr) 4 k 2 + ie J (2vr) 3 J (2tt 



f f d 3 x 2 f d°Xj . 



d 3 3 



(2vr) 3 y (2tt) 



(27T) 



G (W: - k°,q- k)K D {x u q)G Q {W c n - k°,X\ - k)... 
G (W: - k°, x 3 - k)K D (x 3 ,p)G (W: - k°,p- k) 

(-zeY)R c nlj (p,q)} 



— -(Za) 2 [°° dk 

In 7T Jo 



Pi 



P-Q 



d 3 p d 3 q 
(2tt) 3 (2tt) 3 

p 2 + q 2 + 2 



R n i{p)R n i{q) 



(p-q) 2 (k+p 2 + l)(k + q 2 + l) 

1 and n 



From the explicit calculation of ( pj.25| ) on the levels n 
follows: 



(Kladder}-^ = —(Za)" 
71 



| + 81og(2) 



(3.25) 
= 2 it 
(3.26) 
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{Kiadder/ 2 
(K ladder) 



~(Zaf 

7T 



a 



2p 



7T 



(Za) 



497 16, " 

+ — log(2) 

144 3 Sl . 

955 32, . . 

+ — og 2 

432 9 SV ; 



(3.27) 
(3.28) 



These contributions eliminate, as we have announced, the remaining 
a/ir(Za) 2 terms from flOD , 0, O and flpp , O), (p|). 



At the end of this section we present the complete (but for short up to 



order (a/ix) {Za) 7 ) contribution of (K. 
the leading contribution of {K, 



_ . , without the terms cancelled by 

nlj 



rrr / ,. and (Kiadder) 
nlj 



nlj' 



AE ls = (k s . e . + (R ver + K, 





~(Za) 

7T 


+ 


~(Za) 

7T 


+ 


~(Za) 

7T 


+ 


~(Za) 



ladder , . , , 
/ leading I 1 

37 

- + 71og(Za) + 71og(2) 
16 

67T 

9tt 



7T 



1 14 
-/--log 2) 

7T 37T 

983 45 45 

— -_log(Za)- T log(2) 

225tt 4tt 30tt 6V 7 



if. 



-(Za) 

7T 

+ -(Za) 

TT 



+ (k ve 

191 



K, 



ladder 



leading I Is 



256 
5 



+ 
+ 



a 



(Zaf 
(Za) 7 



187T 

2615 



— log (Za) 
16 &V ; 

3tt 1 
T ~ 8tt 



12tt 



log(2) 



12288 
133 



72007T 



509 

256 
3tt 

T 



log (Za) 
53 



128tt~ 



607 
960tt 



log(2) 



(3.29) 



(3.30) 



AE 2p = (k s .e. + (Ever + Ki ad der 



leading I 2p 



16 



-(ZaY 

7T 



43 5 , try ^ 

'768 + 16 l0g(Za) 



+ 



+ 



a 
a 

7i 



{Zaf 
(Za) 7 



9181 329 



36864 
23 



450vr 



768 

7T 1 
- -I 

4 



iog(z«; 



7T 



20vr 



log(2) 



where we have defined []: 

Li 2 (3 - 2^2) + 2Li 2 (l) + log 2 (l + y/2) 



I = 



V2 



(3.31) 



(3.32) 



4 Conclusions 

We have shown how it is possible to perform a perturbative expansion of 
the hydrogenic bound state two-points Green function by choosing an appro- 
priate zeroth-order kernel. We have written the corresponding zeroth-order 
solution and the perturbative expansion for the bound state energy levels. 
Our approach is very transparent and unambiguous in the sense that in this 



7 The presence of the dilogarithms in (3.32) is not surprising, being a standard feature 



of this kind of calculations. For example, such terms can have their origin from integrals 
of the kind 

nl dy 



f(W) = - f ^L\og(y 2 + 2Wy+l), 

Jo y 



obtained after the change of variables p + E p = y and after some rationalizations. It is 
easy to see how 

df(W) 1 VI - W 2 

— r = atari 

dW VI - W 2 W 

IT 2 



/(l) = 2Li 2 (-l) - ( , . 



so that 



6 2 W 

Other terms come from similar integrals giving at the end, after an expansion in Za, the 
results (pj|)-(^2i|). 
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way one knows exactly what he is discarding and what he is keeping, and 
each approximation is just referred to neglecting some Feynman graph. Then 
we have expanded perturbatively the self-energy, and we have calculated the 
exact contribution of the first graph in the expansion, here presented up to 
order (ct/ir) (Za) 7 . The way is still long and hard: the next terms, which 
contribute to the self-energy are the one- Coulomb exchange diagram and the 
sum (from two up to infinity) of all many-Coulomb exchanges. While for 
the one-Coulomb exchange one may think to proceed in the calculation by 
brute force, for the sum of the many- Coulomb exchange graphs there are 
the difficulties to treat the Schwinger integral which compare in the Dirac- 
Coulomb propagator's expansion. These difficulties are essentially the same 
which are present in the other bound-state problems in QED. Therefore, the 
solution of these difficulties is not only important in order to have progress 
in the hydrogenic atoms energy levels calculation but it will make possible, 
in particular, an ex novo calculation in our formalism of the positronium 
energy levels. Along this direction we are going proceed further on. 
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A Appendix 

All graphs of fig. 3 give contributions to the energy level of the same order in 
Za, so that such an expansion of the Dirac-Coulomb propagator gp is wrong 
from the perturbative point of view. In section 2 we have discussed how to 
treat perturbatively in our formalism the Dirac-Coulomb propagator. In this 
appendix we review the alternative method proposed by J. Sapirstein M 
[19| for an application). This method is more artificial than ours which 



sec 



follows as a natural consequence from the formalism. Therefore one expects 
to have some additional analytical problems to implement the following ex- 
pansion for gv into the energy expansion ( |2.7| ). 

The idea is to isolate from ^ a scalar part sd and then to write a perturba- 
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tive expansion for sd- We define sd as: 



g D {W;p,q) = {j°W -p-^ + m)s D (W;p,q) 

f j37 1 



(p — ky 

= (iG l -V cl °)s D . (A.l) 
This means also that so = —ign- 

In the following, in order to make the notation short, we will write all formulas 
as in (|A.1|) , i.e. without the explicit indication of the momenta and their 
integration. Substituting (|A.1| ) into the Dirac-Coulomb equation ( p. 21 ) we 
obtain the equation for sd- 

sd = g {l + ks D ), (A.2) 
k = k c + Sk, (A.3) 



where: 



g (W;p,q) = - — -L^ (2n) 3 6®(p - q), 

W 

k c (W;p,q) = 8ni Za 



(p- q) 2 ' 

5k(p,q) = 2nH(Za) 2 —^— + 4m(Za f ^ ^~ ® 

\p-q\ {p-q} 2 

= Ski + 5k 2 . 

We can consider ( |A.2j ) and ( |A.3| ) as the analogous of (|2.2|), ( |2.3| ) respectively. 
It is important to remark that in ( |A.3|) 5k is really a "perturbation" of k c in 
the sense that all the contributions which arise from 5k are of higher order 
in Za than the contributions which arise from k c . Finally we can write as in 
( |2.6|) the Bethe-Salpeter equation for sj> 



s D = g c (l + 6k s D ), (A.4) 
where g c is given by (see Q2.13Q ): 

1 1 

g c (W;p,q) = g (W;p,q)-— — k c (W;p,q) 



E 2_ W 2 ^ ^^ ! E 2 q -W 2 
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1 + Za 



W 



dp- 



) \/m' 2 -W 2 



^ m 2 _ W 2 J Q ^ p+ ^_ p)2 (E*-WW-W*) 



9o + 9i + 9m- 



A{m 2 -W 2 )(p-q) 2 

(A.5) 



From ( |A.5| ), ( |A.4|) and (|A.1|) we obtain a perturbative expansion for the 
Dirac-Coulomb propagator g D : 



9d 



(iG» 1 - V cl °) 



g c + g c 5k g c + 



+ 9c 5k g c 5k (g c + g c 5k g c + 



(A.6) 



If we want to isolate from this expansion the first two terms of fig. 3, we 
observe that: 

tGu 1 g = (2n) 3 5W(p-q)G , (A.7) 
iG^gl + iG^go 5h g - V cl °g = Go(-ij V c )G , (A.8) 
and hence we rewrite ( |A.6| ) as: 

g D = G + G (-ij°V c )G + iGo l g c M 

+ iGq l g 5k 2 g - V c j°{gt + g c M ) - V c j°g 5k g 

+ (iG^-V c j°) 



(9 C i + 9m) Sk g c + g c 5k {g c x + g c M ) 



+ g c 5k g c 5k (g c + g c 5k g c + ...) 



(A.9) 



We observe that the leading order of the contribution of the remaining graphs 
of fig. 3 arises from iG ~ 1 g M . 
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